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Abstract 

A generic spinning missile with dithering canards is 
used to demonstrate the utility of an overset structured 
grid approach for simulating the aerodynamics of 
rolling airframe missile systems. The approach is used 
to generate a modest aerodynamic database for the 
generic missile. The database is populated with 
solutions to the Euler and Navier-Stokes equations. It is 
used to evaluate grid resolution requirements for 
accurate prediction of instantaneous missile loads and 
the relative aerodynamic significance of angle -of- attack, 
canard pitching sequence, viscous effects, and roll-rate 
effects. A novel analytical method for inter- and 
extrapolation of database results is also given. 

1. INTRODUCTION 

Simulation of the aerodynamics of rolling airframe 
missile systems pose significant challenges for any com- 
putational approach. Applications of practical interest 
are characterized by complex vortical flow and shock 
structures. In addition, the geometry of these missile 
systems can be very complex, involving relative motion 
between missile body and control surfaces. The flows 
are also inherently unsteady. The aim of this paper is to 
demonstrate the utility of Chimera * 1 overset structured 
grid domain decomposition methods in the efficient 
generation of high fidelity aerodynamic simulations for 
this class of problems. 

A generic rolling airframe missile is defined in 
order to demonstrate the advantages and limitations of 
an overset grid approach. The paper provides a techni- 
cal description of the generic missile and the specific 
computational methods here employed, and complete 
discussions of the case conditions and corresponding 
simulation results. The set of simulations considered 
are designed to demonstrate the level of resolution 
required for accurate prediction of surface loads and to 
determine the relative aerodynamic significance of vis- 


cous and missile roll-rate effects. The aerodynamic per- 
formance of the missile as a function of angle-of- attack 
and canard pitching sequence is also considered. A 
novel analytical method for describing these perfor- 
mance characteristics is given. A brief summary, 
acknowledgements, and list of references are provided 
at the end of the paper. 

2. ROLLING AIRFRAME CONFIGURATION 

The generic rolling airframe missile employed in 
the present work is referred to as FM-3. The FM-3 mis- 
sile has a hemispherical nose, cylindrical body, four fins, 
and two canards. Details of the geometric complexity of 
the missile are illustrated in Figure 1. 


a) 




Figure 1 . FM-3 missile geometry, a) Top view of 
entire missile, b) Fins, c) Close-up of canard. 
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The fins are designed to induce missile spin, while 
directional control is actuated via canard dithering. As 
the missile spins, the canard pitch position follows an 
actuator signal with constant pitch rate. The actuator 
signal flip-flops between +/- 1 according to the sign of 
the sum of two sine-waves called the command and 
dither signals. The amplitude of the command signal 
relative to the dither signal is called the command level, 
and reflects the strength of the attempted maneuver. The 
command signal is modulated with the roll-rate. The 
dither signal is modulated with a dither frequency. Fig- 
ure 2 shows the canard pitching algorithm for a com- 
mand level of 100% for the specified roll-rate and dither 
frequency. 

3. COMPUTATIONAL METHODS 
3 . 1 Discretization Method 

The “near-body” and “off-body” domain partition- 
ing method described in References 2 and 3 is used here 
as the basis of discretization of the FM-3 missile. In the 
approach, the near-body portion of a domain is defined 
to include the surface geometry of all bodies being con- 
sidered and the volume of space extending a short dis- 


tance away from the respective surfaces. The 
construction of near-body grids and associated intergrid 
connectivity is a classical Chimera-style decomposition 
of the near-body domain. It is assumed that near-body 
grids provide grid point distributions of sufficient den- 
sity to accurately resolve the flow physics of interest 
(i.e., boundary-layers, vortices, etc.) without the need 
for refinement. This is a reasonable constraint since 
near-body grids are only required to extend a short dis- 
tance away from body surfaces. Figure 3 illustrates the 
surface decomposition of the FM-3 missile and shows 
selected surfaces from the resulting near-body surface 
and volume grids. 

The off-body portion of the domain is defined to 
encompass the near-body domain and extend out to the 
far-field boundaries of the problem. The off-body 
domain is filled with overlapping uniform Cartesian 
grids of variable levels of refinement, as shown in Fig- 
ure 4 for the FM-3 missile. The off-body grid resolution 
amplification factor between successive levels is 2. The 
near-body off-body partitioning approach facilitates grid 
adaptation in response to proximity of body components 
and/or to estimates of solution error within the topologi- 
cally simple off-body grid system. 



Figure 2. Canard pitching algorithm for roll-rate of 8.75 Hz, dither frequency of 35 Hz, and 100% command level. 
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Figure 3. FM-3 surface geometry decomposition and 
near-body grids, a) Surface decomposition of missile 
body. b) Decomposition of canard surfaces and 
selected surfaces from corresponding grids, c) Decom- 
position of fin surfaces and selected surfaces from cor- 
responding grids. 

Note: Surface geometry decomposition and near-body grid genera- 
tion accomplished using OVERGRID utility from CGT 4 . 


Level-3 



x near-body grid 
components 

Figure 4. FM-3 near-body and off-body partitioning 
and selected surfaces. 

3.2 Solution Method 

The set of FM-3 missile simulations presented in 
this paper represent a wide variety of conditions and are 
products of the OVERFLOW-D 3 ’ 5 code. OVERFLOW- 
D is based on version 1.6au of the well known NASA 
OVERFLOW 6 code, but has been significantly 
enhanced to accommodate moving body applications. 
The OVERFLOW-D enhancements represent in-core 
subroutine actuated operations and include the follow- 
ing capabilities. 

i. On-the-fly generation of off-body grid systems. 

ii. MPI 6 enabled scalable parallel computing. 

iii. Automatic load balancing. 

iv. Aerodynamic force and moment computations. 

v. General 6-degrees-of-freedom model. 

vi. Rigid-body relative motion between an arbi- 
trary number of bodies. 

vii. Domain connectivity. 

viii. Solution error estimation. 

ix. Grid adaptation in response to body motion 
and/or estimates of solution error. 

The majority of the FM-3 simulations presented in 
this paper involve relative motion between grid compo- 
nents. The entire missile spins relative to the inertial 
off-body grid system and the canards dither relative to 
the missile body. The pseudo-code below outlines the 
general procedure used in OVERFLOW-D to carry out 
such simulations. Of course, the flow equations are 
solved at every time-step during a simulation. In cases 
that involve relative motion between configuration com- 
ponents, body dynamics and domain connectivity are 
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also addressed at each time-step. In the case of the FM- 
3 missile, body dynamics simply means the computa- 
tion of aerodynamic loads and moving the missile com- 
ponents according to a control-law. Specifically, the 
rotational orientation of the missile is positioned as a 
function of time and the roll-rate. The position of the 
canards relative to the missile body is determined by the 
canard dither algorithm (illustrated in Figure 2). 

For N time-steps 

do every step 

■ Solve flow equations 

■ For Moving Body Problems 

- Body dynamics 

_ - Domain connectivity 

do every m th step 

P i Adaptive Refinement 

- Error estimation 

- Off-body re-partitioning 

- Solution transfer 

|_ - Domain connectivity 

Pseudo-Code. Solution procedure (with adaptive 
refinement capability) for unsteady problems that may 
involve relative motion between component parts. 

Since the missile movement is continuous, the 
relative position of many grid components change every 
time- step. In order for solution information to be cor- 
rectly exchanged between grids during the simulation, 
the domain connectivity solution must also be continu- 
ously updated. This is accomplished automatically by 
OVERFLOW-D. 

The OVERFLOW-D processing rate for static 
geometry viscous flow applications is about 15|Lisec per 
grid-point per time-step (300MHz processor). For mov- 
ing-body problems, the processing rate is somewhat 
problem dependent, but generally falls in the bounds of 
15 to 18 |Llsec per grid-point per time-step. For the FM- 
3 spinning missile cases considered here, the number is 
16.5 jLtsec per grid-point per time-step. OVERFLOW-D 
accommodates problem sizes of more than 2 million 
grid-points per 1 gigabyte of memory. Maximum paral- 
lel efficiency (percentage in high 90’ s) is realized when 
the fewest number of processors that can accommodate 
a given problem in core memory are selected. OVER- 
FLOW-D can efficiently (i.e., over 70%) make use of 
larger numbers of processors for a fixed problem size 


when each processor assumes the load of at least 250 
thousand points. Load balancing is an automatic func- 
tion of OVERFLOW-D. 

As indicated in the pseudo-code above, OVER- 
FLOW-D accommodates solution adaptation based on 
the position of near-body grid components and/or in 
response to estimates of solution error. The off-body 
grid management scheme allocates level- 1 (finest) reso- 
lution grids to accommodate significant motion of body 
components or flow features before the next adapt cycle. 
Accordingly, adapt cycles are only required periodi- 
cally; every 25 to 50 time-steps in a typical unsteady 
simulation. 

In the FM-3 missile cases considered here, all flow 
features that are likely to have any significant affect on 
the surface forces and moments are confined to the vol- 
ume of space within a missile diameter of the body 
itself. These include canard vortices, boundary layer, 
and key portions of the shock systems. Accordingly, 
OVERFLOW-D input is used to allocate level- 1 resolu- 
tion capacity to a distance of 1.5 diameters from the 
missile surface, rather than enable adaptation in 
response to solution error. A slight savings in computa- 
tional overhead is thereby gained for the present cases. 

4. SIMULATION RESULTS 

A set of FM-3 missile simulations is carried out to 
demonstrate the level of resolution required for accurate 
prediction of surface forces and moments and to deter- 
mine the relative aerodynamic significance of viscous 
effects, missile roll-rates (O r ), canard command levels 
(c), and free-stream angles-of-attack (a). A total of 31 
FM-3 simulations form the basis of the material pre- 
sented in this paper concerning these issues. The flight 
conditions for the cases are indicated in Table 1. The 
parameters varied to obtain the complete simulation set 
are free-stream angle-of-attack, canard pitch command 
level, and missile roll-rate. The other simulation param- 
eters indicated in the table are held fixed and are com- 
mon to all cases considered. 


Table 1 . Simulation Parameters 


Moo 

Re 

Mach number 
Reynolds number 
Dither-frequency 

1.6 

50 x 10 6 
35 Hz 

a 

Angle-of-attack 

0°, 2°, 3°, 4°, 8°, 12°, 15° 

c 

Command level 

0%, 100%, 200% 

O r 

Roll-rate 

0 Hz, 8.75 Hz 
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The general characteristics of the FM-3 flow field 
are illustrated in Figure 5. Vortices are shed from the 
inboard and outboard tips of the canard and convect 
down the length of the missile interacting with the vis- 
cous boundary layer. Away from the influence of the 
boundary layer, the outboard canard vortices twist 
approximately 8° around the spinning missile in one 
body-length of travel. As can be seen in the figure, dis- 
ruptions to the boundary layer by the inboard canard 
vortices are dragged through nearly 45° of roll over the 
same interval. The shock structure is indicated in Figure 
5b. The boundary layer growth on the upper surface of 
the missile is also visible in the figure. A positive angle- 
of-attack and vortex/boundary layer interaction combine 
to exaggerate the boundary layer thickness down the 
stream- wise axis of the missile. 



Inboard canard 
vortex 


Outboard canard 
vortex 


Vortex/Boundary layer 
interaction 



b) 


Figure 5. Aerodynamics of the FM-3 spinning missile 
with dithering canards, a) Vortex structure, b) Shock 
structure. £2 r = 8.75 Hz, c = 0%, a = 3°. 


4.1 Resolution Requirements 

A grid refinement study is used to determine the 
level of spatial resolution needed to accurately predict 
the integrated FM-3 surface loads. The significance of 
grid resolution is evaluated here by comparing viscous 
solutions for the spinning missile case defined in Table 1 
with the variable parameters c, a, and O r fixed at 0%, 


3°, and 8.75 Hz, respectively. A very high resolution 
grid is used to define the baseline solution. Medium and 
coarse solutions are obtained on grids derived from the 
baseline grid with successively lower levels of spatial 
resolution. The qualitative effect of coarsening on the 
surface geometry is shown in Figure 6. A very high 
fidelity temporal resolution (viz., 12,000 time-steps per 
missile revolution) is uniformly employed in all of the 
viscous simulations. 



Figure 6. Fine, medium, and coarse grid representa- 
tions of FM-3 canards and fins. 


The baseline grid for this case (finest resolution) is 
comprised of 41 million grid points and is referred to 
hereafter as the VI (i.e., “Viscous-1”) grid. Isolated sur- 
faces from the VI grid are shown in Figures 3 and 4. 
The flow characteristics illustrated in Figure 5 are from 
a simulation using the V 1 grid. All grid lengths referred 
to in the following discussion are normalized by the 
missile body length. Viscous spacing normal to the 
body surfaces is 2.5 x 10 in the VI grid. This corre- 
sponds to a y + of 1 for a Reynolds number of 10 million. 
This spacing is maintained uniformly across the first 6 
cells in the viscous direction and then expanded with a 
geometric stretching ratio of 1.2 to a distance of approx- 
imately 0.015. The maximum spacing used in the near- 
body grids is approximately equal to the level- 1 off- 
body grid spacing which is 0.0013, or approximately 
0.1% of the body length. 
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The VI grid is the basis of the medium (V2) and 
coarse (V3) grids. The V2 grid is obtained by deleting 
approximately every other point from the V 1 grid in all 
three spatial dimensions and results in a grid with just 
over 8 million points. Similarly, the V3 grid is obtained 
from the V2 grid by deleting approximately every other 
point from the V2 grid in all three spatial dimensions 
and results in a grid with just over 2 million points. The 
foregoing is true subject to the following qualifications. 

i. Some surface grids require redistribution and / 
or addition of grid-points to preserve geometric 
features such as sharp corners through the two 
subsequent eliminations of every other point. 

ii. Smoothing is applied to geometric features that 
are not adequately resolved by the coarser 
grids. For the V2 grid, smoothing is applied to 
the canard cut-out and the missile groove. For 
the V3 grid, the canard cut-out, the canard 
shaft and the missile groove are removed com- 
pletely. 

iii. The grid spacing in the surface normal direc- 
tion for the V2 grid corresponds approximately 
to every other point for the V 1 grid, doubling 
the initial spacing from the wall. 

iv. The grid spacing in the surface normal direc- 
tion for the V3 grid starts at the surface with 
the V2 spacing doubled. The stretching ratio 
thereafter is approximately the same as in the 
V2 grid. 

Consider now the computed load histories obtained 
from simulations using the VI, V2, and V3 grids. Fig- 
ure 7 shows the computed normal force (Cz) history for 
the three different resolution capacities. The canard 
pitch angle history is also indicated. The VI and V2 
results are in good agreement, except at maximum 
canard deflection. At high canard lift, the strong vorti- 
ces shed from the canards modify the pressure distribu- 
tion on the fuselage and the tail fins. Still, the roll- 
averaged normal forces from V 1 and V2 shown in Table 
2 differ by less than 0.3%, indicating near grid conver- 
gence for this quantity. The V3 result differs signifi- 
cantly from V 1 and V2. 

The data represented in Figures 8 and 9 observe the 
same form as that used in Figure 7. However, Figures 8 
and 9 display side force (Cy) and axial force (Cx) histo- 
ries, respectively. The side forces exhibit the same 
effect as for the normal forces at maximum canard 
deflection. The roll-averaged side-forces are close to 
zero, with a difference between VI and V2 of less than 
0.02% of the normal force, indicating grid convergence 
for this quantity also. In contrast, the axial forces (see 


Table 2. Roll-averaged force and moment coefs.* 


Coefficient 

Fine (VI) 

Med. (V2) 

Coarse (V3) 

Cx (axial) 

1.17 

1.12 

1.07 

Cy (side) 

-7.56e-03 

2.76e-03 

-1.21e-03 

Cz (normal) 

0.461 

0.462 

0.538 

Cmx (roll) 

-1.19e-03 

-1.08e-03 

-0.94e-03 

Cmy (pitch) 

-6.79e-03 

-3.11e-03 

-3.45e-02 

Cmz (yaw) 

2.80e-03 

3.52e-03 

-9.55e-04 


Moments are about the missile center of gravity 


Figure 9), indicate a systematic shift higher for finer res- 
olution. The difference between VI and V2 axial 
forces is approximately the same as between the corre- 
sponding V2 and V3 results. Clearly, grid convergence 
is not apparent in the computed axial force data. 

Figures 10 and 11 break down the axial forces into 
pressure and viscous components. The V2 and V3 solu- 
tions are almost identical for the pressure component; 
and VI has a systematic shift to a higher value. The 
contribution to this shift comes mainly from the aft part 
of the missile. This region of the flow has complicated 
interactions between the expansion waves around the 
boat-tail, shocks around the tail fins, and the boundary 
layer. The VI and V2 viscous components to the axial 
force differ by approximately 2% of the total axial force, 
slightly less than the difference between the V2 and V3 
results. The total roll-averaged axial forces are shown 
in Table 2. The medium and coarse solution axial forces 
are 4% and 8% lower than for the corresponding fine 
solution. The reason grid convergence in axial force is 
not demonstrable via the current set of solutions is not 
clear. It may be that while the V 1 grid has sufficient res- 
olution in the boundary layer, the V2 and V3 viscous 
spacing (double and quadruple that of VI) is not suffi- 
cient. 

Figure 12 shows the pitching moment (Cmy) about 
the center of gravity for the VI, V2, and V3 simulations. 
The missile center of gravity is located approximately at 
the missile midpoint. As is the case for normal force, 
the differences in pitching moment are largest at maxi- 
mum canard deflection. Still, the overall agreement is 
very good. The difference between the VI and V2 roll- 
averaged pitching moments are approximately 1% of the 
maximum pitching moment during a revolution. The 
roll-averaged pitching moments are shown in Table 2. 
The percentage of maximum pitching moment is used 
here as a measure of grid convergence since the roll- 
averaged moments are all nearly zero. 

Figure 13 shows the yawing moment (Cmz) about 
the center of gravity. The agreement is very good, with 
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Figure 7. Grid effects in the NORMAL FORCE coefficient. 



Figure 8. Grid effects in the SIDE FORCE coefficient. 
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Figure 9. Grid effects in the AXIAL FORCE coefficient. 



Figure 10. Grid effects in the AXIAL FORCE coefficient (PRESSURE component). 
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Figure 11. Grid effects in the AXIAL FORCE coefficient (VISCOUS component). 
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Figure 12. Grid effects in the PITCHING MOMENT coefficient. 
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Figure 13. Grid effects in the YAWING MOMENT coefficient. 


some minor differences at maximum canard deflection. 
The difference between the VI and V2 roll-averaged 
yawing moments are approximately 0.2% of the maxi- 
mum pitching moment during a revolution. The roll- 
averaged yawing moments are shown in Table 2. 

The rolling moment (Cmx) is shown in Figure 14. 
As expected, the moments are small. The difference 
between the VI and V2 roll-averaged rolling moments 
are approximately 0.03% of the maximum pitching 
moment during a revolution. A steady-state freely spin- 
ning missile should have an averaged rolling moment of 
zero. The overall pattern here is in fair agreement, but 
the details differ significantly. The rolling moment is 
sensitive to flow details around the tail fins, with com- 
plex interactions between the boundary layer, the canard 
vortices and the shocks, which again depend on the 
accurate prediction of the flow along the entire missile. 
Grid convergence for so small a quantity as the rolling 
moment is beyond the capacity of the present set of 
grids. 

The effect of grid resolution on the prediction of the 
aerodynamic details of the flow is illustrated in Figure 
15. Figures 15a, 15b, and 15c provide a comparative 
view of the canard vortices and missile boundary layer 
interaction via plots of helicity density (i.e., dot product 
of the velocity and vorticity vectors) at several stations 
along the length of the missile (for a roll-angle of 268.5° 


at maximum negative canard deflection). The position 
of the vortices are in good agreement for the V 1 and V2 
solutions. However, the vortex strength is weaker and 
interactions between the inboard canard vortex and vis- 
cous boundary layer of the missile body are less appar- 
ent in the V2 solution. The V3 solution differs 
significantly from the VI and V2 solutions in vortex 
position, strength, and vortex/boundary layer interac- 
tion. 

The grid refinement results suggest very good over- 
all agreement between the VI and V2 solutions, though 
some differences are apparent. Still, the V2 grid offers a 
good compromise between solution accuracy and solu- 
tion throughput for computations designed to predict 
aerodynamic forces and moments. As noted above, grid 
convergence of all forces and moments (except axial 
force and rolling moment) are obtained. 

4.2 Viscous Effects 

A comparative evaluation of very high resolution 
Navier- Stokes and Euler simulations is used to deter- 
mine the relative significance of viscous effects opera- 
tive in the range of flight conditions considered for the 
FM-3 missile. The case conditions defined in Table 1 
are taken as representative of these flight conditions. A 
static FM-3 case is first considered where the missile 
roll-rate is zero and the canards are fixed in neutral posi- 
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Figure 14. Grid effects in the ROLLING MOMENT coefficient. 



b) V2 (medium) solution 


c) V3 (coarse) solution 


a) V 1 (fine) solution 


Figure 15. Fine, medium, and coarse grid solutions for the FM-3 missile (snap-shot at a roll-angle of 268.5°). 
M x = 1.6, Re = 50xl0 7 , a = 3 °,c = 0%, O r = 8.75 Hz, O d = 35 Hz. 
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tion. The fundamental differences between static vis- 
cous and inviscid solutions are apparent and are 
germane to spinning cases. A spinning FM-3 missile 
with dithering canards case is then considered in detail 
for a command level of 0%, free-stream angle-of-attack 
of 3°, and missile roll-rate of 8.75 Hz. 

The VI solution considered in Section 4.1 is used 
as the basis of comparison for a correspondingly high 
resolution inviscid simulation. The high resolution 
inviscid, or “Euler grid,” is referred to hereafter as the 
El grid system. The El surface grids and off-body vol- 
ume grids are identical with the corresponding compo- 
nents of the VI grid system. The only differences 
between the El and VI grids is the surface normal dis- 
tribution of points in the respective near-body grid com- 
ponents. The El surface normal wall spacing is 20 
times that of the VI grid. 33 million points are used to 
define the complete El grid compared to the 41 million 
points used in V 1 . 

The computational savings available by assuming 
inviscid flow are significant. In the present high resolu- 
tion cases, 20% fewer grid points are used in the El sys- 
tem than in the VI system. Due to the larger surface 
normal wall-spacing, larger stable time-steps are also 
possible - a At increase of 5 times is used in the present 
El simulations, allowing for nearly 2,500 steps per mis- 
sile revolution. In the present simulations, the combined 
effects of fewer grid points, fewer floating-point opera- 
tions required per grid point, and larger time-steps result 
in an order of magnitude savings in computational 
expense. 

4.2.1 Static Geometry FM-3 Case 

Consider the qualitative differences between vis- 
cous and inviscid FM-3 missile solutions for zero roll- 
rate and neutral canard positioning. These differences 
characterize some of the trends that are apparent for the 
spinning missile and canard dithering conditions of 
interest here. The vortex structure of the viscous and 
inviscid non-spinning FM-3 missile solutions are shown 
in Figure 16. The corresponding shock structure and 
surface pressure distributions for the two cases are 
shown in Figures 17 and 18, respectively. 

The position and strength of the outboard canard 
vortices are essentially identical in both solutions. 
However, the inboard canard vortices differ significantly 
in both strength and position. Opposed by the viscous 
action of the boundary layer, the inboard canard vortices 
of the viscous solution are weaker and traverse a differ- 
ent path than their inviscid counterparts. The inboard 


canard vortices (viscous case) also pull the boundary 
layer off the missile surface near the tail. 

At the missile nose, the shock structure is similar 
for the viscous and inviscid solutions. The main fea- 
tures are a normal shock in front of the nose, an oblique 
shock above and below the nose, expansion around the 
two sharp corners between nose and fuselage, and the 
shocks in front of the canards. The local influence of 
the canard on the fuselage pressure distribution is 
smeared somewhat in the viscous case, due to the 
boundary layer. The expansion and compression over 
the groove is similarly weaker in the viscous case due to 
the boundary layer. At the tail, the interaction between 
the boundary layer and the expansion- waves and shock- 
waves influences the pressure distribution significantly. 

4.2.2 Dynamic Geometry FM-3 Case 

Consider viscous and inviscid solutions for the 
spinning FM-3 missile (£2 r = 8.75 Hz and = 35 Hz). 
The corresponding roll-averaged force and moment 
coefficients are given in Table 3. Figure 19 shows the 
normal force histories over one revolution. The viscous 
and inviscid computations agree well, except at maxi- 
mum canard deflection, where there are instantaneous 
differences of up to 20%. The influence of the deflected 
canards on overall forces appear more pronounced in the 
inviscid case, though this is an indirect effect. Figures 
20, 21, and 22 decompose the normal force into contri- 
butions from the canards, fuselage, and tail section, 
respectively. There is almost no difference between the 
viscous and inviscid canard normal force histories 
shown in Figure 20. The over prediction of normal 
forces at maximum canard angles in the inviscid solu- 
tion are realized from pressure distribution differences 
along the fuselage and tail sections. The vortex systems 
shed from the canards interact with the fuselage and tail 
fins differently in the viscous and inviscid solutions. 
The inviscid inboard canard vortices are too energetic 
and positioned incorrectly, leading to the differences 
seen in Figures 21 and 22. The same effects are true for 
the side forces shown in Figure 23. 

Table 3. Roll-averaged force and moment coefs.* 


Coefficient 

Viscous (VI) 

Inviscid (El) 

Cx (axial) 

1.17 

1.10 

Cy (side) 

-7.56e-03 

-19.6e-03 

Cz (normal) 

0.461 

0.449 

Cmx (roll) 

-1.19e-03 

-1.65e-03 

Cmy (pitch) 

- 6.79e-03 

-1.06e-03 

Cmz (yaw) 

2.80e-03 

9.17e-03 


Moments are about the missile center of gravity 
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Figure 16. Vortical structure for static geometry FM-3 missile. M ^ = 1.6, Re = 50xl0 7 , a = 3°. 


viscous 


inviscid 


Figure 17. Shock structure for static geometry FM-3 missile. M ^ = 1.6, Re = 50xl0 7 , a = 3°. 
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a) viscous 



c) d) 


Figure 18. Surface Cp distributions for static geometry FM-3 missile. = 1.6, Re = 50xl0 7 , a = 3°. 



Figure 19. Viscous effects in the NORMAL FORCE coefficient. 
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Figure 20. Viscous effects in the NORMAL FORCE coefficient (CANARD COMPONENT ONLY). 



Figure 21. Viscous effects in the NORMAL FORCE coefficient (FUSELAGE COMPONENT ONLY). 
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Figure 22. Viscous effects in the NORMAL FORCE coefficient (TAILFIN COMPONENTS OMA 



Figure 23. Viscous effects in the SIDE FORCE coefficient. 
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The reason that the viscous and inviscid solutions 
agree very well everywhere except at maximum canard 
deflection is evident in the vortical structures shown in 
Figure 24. Figure 24a shows cutting planes with helic- 
ity density contours at several stations along the length 
of the missile for a roll-angle of 61.5°. At this roll- 
angle, the canard is in transition between the maximum 
positive and negative canard pitch angles. Vortices trail 
each canard in counter-rotating pairs. Whenever the 
canards snap from positive to negative (or negative back 
to positive) pitch angle, the rotational sense of the vorti- 
ces also reverse direction. 

This effect is clearly evident in the instantaneous 
snap-shot of the vortical field shown in Figure 24a. 
From the canard downstream to about the missile mid- 
point, the vortices become weaker and finally disappear. 
From the midpoint back to the tail section, the vortices 
re-appear rotating in the opposite sense and growing 
slightly in strength. The history of one transition cycle 
is captured in this image. Vortex interaction with the 
boundary layer is minimal during canard transition. 
Accordingly, the viscous and inviscid solutions are in 
very good agreement for all roll-angles where the 
canards are in transition. 

Figure 24b shows cutting planes with helicity den- 
sity contours at several stations along the length of the 


missile for a roll-angle of 84°. The instantaneous solu- 
tion indicated in the figure corresponds to the comple- 
tion of nearly 20 degrees of missile role, or 2.5 missile 
body-lengths of travel, with the canards positioned at 
maximum negative deflection. The fully developed 
canard vortices prevail well past the missile. The vorti- 
ces interact with each other and with the boundary layer. 
The inboard canard vortices are weakened by the bound- 
ary layer, reducing the force peaks, as seen in the normal 
and side force histories of Figures 19 and 23. The 
inboard canard vortex/boundary layer interactions are 
entirely unaccounted for in an inviscid simulation. 

The roll-averaged axial force for the viscous solu- 
tion is about 8% higher than the corresponding inviscid 
solution (see Table 3). Figure 25 shows the correspond- 
ing axial force histories over one revolution. Approxi- 
mately 60% of the difference noted between the viscous 
and inviscid axial force is due to viscous wall stresses. 
The remaining 40% of the difference is due to spatial 
variations in the surface pressure distributions. The 
pressure contribution to axial force in the viscous solu- 
tion is compared separately in Figure 26 with the invis- 
cid axial force. Boundary layer/shock interactions in 
the viscous solution lead to differences in overall shock 
structure and surface pressure distributions, accounting 
for these differences. 
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b) Maximum negative canard deflection, 0 = 84° 


Figure 24. Vortical structure for vsicous FM-3 spinning missile case. M ^ = 1.6, Re = 50xl0 7 , a = 3°, c = 0%, 
£l r = 8.75 Hz, £2 d = 35 Hz. 
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Figure 25. Viscous effects in the AXIAL FORCE coefficient. 
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Figure 26. Viscous effects in the AXIAL FORCE coefficient. 
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Figure 27 shows the pitching moment about the 
missile center of gravity. The viscous solution has 
higher peaks. The stronger vortices in the inviscid solu- 
tion interact with the fuselage and tail fins, adding a 
moment that counteracts the moment from the canards. 
A similar effect can be seen for the yawing moment in 
Figure 28. 

Based on the high resolution viscous (VI) and 
inviscid (El) simulations, a few statements regarding 
the relative significance of viscous effects on FM-3 mis- 
sile performance are justified. 

i. Instantaneous side- and normal forces differ up 
to 20%. This is due to differences in interaction 
between the vortex system shed from the 
canards and the fuselage boundary layer and 
tail fins. 

ii. The axial forces differ by about 8%. The vis- 
cous wall stresses accounts for 60% of this dif- 
ference. 

iii. Euler computations provide valuable informa- 
tion, but viscous effects should ideally be 
included in detailed studies. 

4.3 Roll-Rate Effects 

An important objective in considering the utility of 
computational methods (and physical experiments) for 
rolling airframe missile systems is the accurate predic- 
tion of roll-averaged forces and moments. A large 
parameter space can easily exhaust computational (or 
experimental) resources, making it important to elimi- 
nate parameters of secondary importance. The spinning 
missile cases defined in Table 1 have a roll-rate of 8.75 
Hz, which allows for only 8° of missile roll per body- 
length of travel. It seems plausible that for a given 
pitching sequence and moderate angular velocity, the 
roll-averaged forces may not be sensitive to roll-rate. 
The relative significance of roll-rate effects in this 
regime are determined here by evaluating differences 
between solutions having roll-rates of zero and 8.75 Hz. 
The baseline flight conditions defined in Table 1 and the 
V2 grid system are used in all cases. The canard com- 
mand level is 0% and the free-stream angle-of-attack is 
3° 

The 8.75 Hz spinning missile case is discussed at 
length in Section 4.1 and results using the V2 grid sys- 
tem are given in Table 2 and Figures 7 through 15. The 
“zero” roll-rate cases considered in the present section 
for comparative purposes are generated in a discrete 
quasi-static fashion for the same roll-angle and canard 
dither cycle. The quasi- static results for the zero roll- 


rate condition can be obtained in one of several ways. 
Two obvious ways are as follows. 

i. Static geometry . Start with the desired missile 
roll and canard pitch angles fixed in free- 
stream conditions. Drive the corresponding 
solution to convergence in steady-state mode. 
Switch to time-accurate mode (static geome- 
try) and continue the solution to determine the 
temporal state of the flow. 

ii. Frozen geometry. Freeze the geometry of the 
8.75 Hz spinning case at the desired missile 
roll and canard pitch angles. Initialize the solu- 
tion with the corresponding instantaneous spin- 
ning result obtained previously (see Section 
4.1). Hold the roll and canard pitch angles 
constant and run the simulation time-accu- 
rately until the false-transient dies and to deter- 
mine the temporal state of the flow. 

Figure 29 shows the time-history of the normal 
force that results from the “static geometry” and “frozen 
geometry” approaches corresponding to a roll-angle of 
150° and -6.95° canard pitch angle. The normal force 
coefficient is plotted versus “relative” time-step, 
because both steady-state and time-accurate integration 
are used during the convergence histories. For the static 
geometry case, the solution is run to convergence in 
steady-state mode. At relative time-step zero, the time- 
integration scheme is switched to time-accurate mode 
and run further as shown. In all, more than 15,000 time- 
steps are needed to obtain a converged solution with this 
approach. The number of time-steps (or iterations) can 
be reduced with alternative integration schemes (e.g., 
dual time-stepping, multigrid, etc.), but the effort is still 
a significant fraction of the cost of computing a full rev- 
olution for a dynamic case. 

In the frozen geometry result shown in Figure 29, 
the instantaneous solution from the spinning missile 
case is used for initial conditions. At relative time-step 
zero, the missile orientation is frozen to the identical 
state used with the static geometry approach and the 
solution is commenced time- accurately. Cz quickly set- 
tles (in less than 500 steps) to a value that is identical to 
the Cz obtained from using the static geometry 
approach. Time histories for the other forces and 
moments show similar behavior. A comparison of the 
two methods can be summarized as follows: 

i. The two methods provide the same results for 
forces and moments. 

ii. A restart from the dynamic flow-field provides 
a steady- state solution at least 20 times faster 
than starting from scratch for this particular 
roll-rate. 


American Institute of Aeronautics and Astronautics 

19 




Figure 27. Viscous effects in the PITCHING MOMENT coefficient. 



Figure 28. Viscous effects in the YAWING MOMENT coefficient. 
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Figure 29. NORMAL FORCE coefficient. Comparison of “static geometry” and “frozen geometry” approaches to 
generating quasi-static solutions for rolling airframe missiles. Static geometry converged in steady-state mode until 
step 0, then run time- accurately. Frozen geometry case run time-accurately from spinning missile initial condition. 

ill. Even if the problem turns out to be quasi- vorticity dynamics become significant. Since the VI, 

steady, it is more efficient to run through a V2, and V3 grid resolution study does not demonstrate 

dynamic computation rather than computing a convergence for the rolling moment, an analysis of roll- 

series of static solutions. rate effects on rolling moment is not attempted here. 

The difference between the frozen geometry initial In a study of the significance of roll-rate effects on 

solution and the converged static geometry solution the FM-3 missile aerodynamics, six data points can only 
directly shows the point-wise relative effect of roll-rate. give an indication of the complete phenomenon. Still, it 

Figure 30 shows a comparison of 6 static solutions seems plausible that moderate changes of roll-rate in the 

(obtained via the frozen geometry method) to the time- range of 0 - 10 Hz should not affect roll-averaged forces 

accurate solution for the 8.75 Hz spinning case. The significantly, 

corresponding static and dynamic results are in good 

agreement. For roll angles 30°, 90°, 210° and 270°, the 4.4 Computational Expense 
dynamic results lag the static results by a few degrees of 

roll angle. For the roll angles 150° and 330°, the canard The FM-3 missile cases converge to a periodically 

is in the middle of a canard "snap" between the two repeating solution in approximately 460° of roll. The 

extreme deflections. This causes a dynamic inflow computational expense of each V2 time-accurate solu- 

effect, with both a delay and larger changes of forces on tion is approximately 438 hours of CPU time on a 300 

the canards than what is sustained in a static configura- MHz processor (viz., SGI Origin 2000) per revolution, 

tion of the missile. At 150° roll angle, the force over- or 560 CPU hours per case. Most of the V2 results pre- 
shoot outweighs the lag effect, and shifts the dynamic sented in this paper are the result of runs using 16 pro- 
solution ahead of the static solution. At a roll angle cessors, with a processing rate of 35 hours per case and 

330°, the effects seem to cancel. approximately 95% parallel efficiency. 

The other forces and moments (except rolling The OVERFLOW-D performance rate realized for 

moment) show similar behavior. The rolling moments the V2 solutions justifies the contemplation of aerody- 

differ significantly. This is not surprising, since the roll- namic database population with Navier-Stokes solu- 

ing moment is generated by the tail fins for which the tions. Maximum parallel efficiency per case is realized 
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Figure 30. 



when the minimum number of processors that will 
accommodate the case in core memory is used. Four 
processors are sufficient for the V2 solutions (assuming 
1 gigabyte of memory per processor). Simultaneous 
computation of multiple cases leads to perfect parallel- 
ism in the generation of solution sets. For example, a 
256 node SGI Origin 2000 can execute 64 FM-3 cases 
simultaneously yielding a throughput of approximately 
330 cases per month. A database populated by conven- 
tional means (i.e., physical experiments) does not 
require many more points than this. For a rolling air- 
frame, such as the FM-3, approximately 600 spinning 
missile data points is sufficient for a parameter space 
that includes Mach number, angle-of-attack, command 
level, and roll-rate. 

5. AERODYNAMIC PERFORMANCE 
5 . 1 Simulation Results 

Given an understanding of the grid resolution 
requirements and the relative significance of viscous and 
roll-rate effects, a matrix of cases is defined to evaluate 
the aerodynamic performance of the FM-3 spinning 
missile subject to variations in free- stream angle-of- 
attack (a) and canard command level (c). The baseline 
conditions for the case matrix are defined by the fixed 
parameters of Table 1 (M m , Re, and O d ) and with the 


roll-rate (O r ) set to 8.75 Hz. The matrix consists of 21 
unique cases defined by varying a and c over the ranges 
indicated in Table 1 (viz., a = 0°, 2°, 3°, 4°, 8°, 12°, 15° 
and c = 0%, 100%, 200%). The pitching sequences that 
result from the three command levels are shown in Fig- 
ure 31. 

The results given in Section 4 demonstrate that 
although inviscid computations can provide valuable 
information about the FM-3 spinning missile, there are 
significant viscous effects (e.g., fuselage boundary layer 
damping of the inboard canard vortex). The results also 
indicate that the medium resolution viscous grid (V2) 
yields comparable results to the 41 million point base- 
line viscous grid (VI). Accordingly, the case matrix 
considered here is populated entirely with viscous solu- 
tions using the V2 grid system. The results of Section 4 
also suggest that roll-rate effects are not important for 
the parameter space of interest here. However, the cost 
of generating 21 time-accurate spinning missile solu- 
tions to populate the case matrix is far less expensive 
than generating a comparable data base composed of 
quasi-static solutions. Time-accurate simulation data is 
therefore used. 

The normal forces for the 7 cases with command 
level 0% are shown in Figure 32. The corresponding 
results for the 100% and 200% command level are 
shown in Figures 33 and 34, respectively. Results for 
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Figure 31. Canard pitching sequences for 0%, 100%, and 200% command level. 



Figure 32. NORMAL FORCE coefficient for different a at 0% command level. 
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Figure 33. NORMAL FORCE coefficient for different a at 100% command level. 



Figure 34. NORMAL FORCE coefficient for different a at 200% command level. 
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the 2° and 4° angle-of-attack fall smoothly between 
neighboring data points, and are omitted in the figures 
for clarity. The effect of increased a at each command 
level is a shift of the normal force curves upwards, with 
a limited change in variation with roll-angle. The direct 
contribution of the canards to instantaneous normal 
force dominate all other sources. However, the fuselage 
and tail section contributions are not insignificant (recall 
Figures 20 through 22). Still, the variation in fuselage 
and tail-section contributions to normal force are a sec- 
ondary effect of canard position and resulting vortex/ 
body interaction. 

The side forces for the 7 cases at 0% command 
level are shown in Figure 35. Again, results for a = 2° 
and 4° are omitted for clarity. The effect of a here is 
less clear than for the normal force. An increase of a 
increases the variation of side force. An interesting 
result apparent in Figure 35 are non-zero side forces at 
roll-angles of 0° and 180° (canards horizontal). This is 
due to the induced circulation from the missile spin 
interacting with the vertical component of the upstream 
velocity vector. This is similar to the Magnus effect on 
a spinning golf-ball. 

The roll-averaged normal forces are shown in Fig- 
ure 36 and, consistent with the trend in Figures 32 


through 34, increase with command level. The roll- 
averaged side forces are shown in Figure 37. The side- 
forces increase with command level. Further, at each 
command level, the side force goes from positive to neg- 
ative with increasing angle-of-attack. The negative 
side-force at high a is caused partly by the Magnus 
effect as mentioned earlier. For a rotating cylinder, the 
Kutta-Joukowski theorem can be evaluated as follows: 

C y = -4sinagf (1) 

where a is the angle-of-attack, £1 the cylinder spin-rate, 
L the cylinder length, and U ^ the speed of the free- 
stream flow. The result is shown in Figure 37 for com- 
parison. 

The roll-averaged pitching moments about the cen- 
ter of gravity are shown in Figure 38. The increase in 
pitching moment with command level is consistent with 
the increase in normal force for the canard. A change in 
command level to maneuver the missile affects the 
pitching moment in a similar manner for different a. 

The roll-averaged yawing moments about the center 
of gravity are shown in Figure 39. The negative contri- 
bution from an increase of command level is consistent 
with the increase in canard side-force. The effect of 



Figure 35. SIDE FORCE coefficient for different a at 0% command level. 
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Figure 36. Roll-Averaged NORMAL FORCE coeffi- 
cients for the Case Matrix. 



Figure 38. Roll- Averaged PITCHING MOMENT coef- 
ficients for the Case Matrix. 


induced circulation seen on the side-force hardly modi- 
fies the yawing moment about the center of gravity, 
since it is relatively evenly distributed along the missile 
length. 

5.2 Interpolation of Results 

The roll-averaged inertial-frame forces and 
moments are the object of analysis in preceding sections 
of this paper. Indeed, the unsteadiness seen in instanta- 
neous results is primarily caused by motion of the 
canards. However, comparison between canard motion 
and resulting forces and moments is complicated by the 
rotating framework used for presentation of the pitching 
sequence. It is not obvious, for example, that the com- 
puted force in Figure 30 corresponds to the pitching 



Figure 37. Roll- Averaged SIDE FORCE coefficients for 
the Case Matrix. 



Figure 39. Roll- Averaged YAWING MOMENT coeffi- 
cients for the Case Matrix. 


sequence in the same figure. Looking at the three pitch- 
ing sequences in Figure 31, it is not immediately clear 
what to expect of changes in roll-averaged forces when 
the command level is changed, and if these changes can 
be modeled well by linear interpolation, for example. 

A method of classifying a canard pitching sequence 
is proposed here to aid in the understanding of com- 
puted results, and to provide a tool for interpolation of 
roll-averaged forces and moments between different 
command levels. The approach taken is to derive a mea- 
sure of the aerodynamic force on the canards as a func- 
tion of canard orientation, which is basically the same 
approach as using angle-of-attack to evaluate lift on a 
2D airfoil. Here, the force (vector) is averaged over one 
missile revolution to represent the mean force vector on 
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the canards. As discussed earlier, the canard forces are 
not the only contributors to changes in overall missile 
forces, but should be the best single-parameter approxi- 
mation. 

Two cartesian right-handed co-ordinate systems are 
used in the analysis. The x-y-z system is in the inertial 
frame, with the x-axis pointing nose to tail, y-axis point- 
ing to the side, and the z-axis pointing up (see Figure 
40). The r-s-t system is initially oriented as the x-y-z 
system, corresponding to horizontal canards with zero 
pitch. However, the r-s-t system rotates with the 
canards. The movement of the r-s-t basis vectors are 
evaluated as follows: 

i. The pitch axis e s is rotated around the inertial 
frame x-axis with the roll-angle. 

ii. The ~e r basis vector is rotated around the new 
e s with the canard pitch angle. 

iii. The e t basis vector is computed as the cross- 
product e r xe s 



e t 


Figure 40. r-s-t rotating reference frame basis vecors 
relative to inertial frame x-y-z coordinate axes. 


The canard angle-of-attack vector, a c , is proposed 
in Equation 2 as a single parameter approximation for 
changes in missile loads due to changes in canard com- 
mand level. 


2 2 
W, + u f 


sin ( 
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u is the upstream velocity vector at the missile angle-of- 
attack, a. It has components u r , u s , and u t in the rotat- 
ing system. 

The three terms in Equation 2 are: 

1 A direction vector normal to the canard. This 
approximates the direction of the force acting 
on the canard. 

2 A fraction of the dynamic pressure, varying 
between 0 and 1 . This reduces the force on the 
canard according to the local "side-slip" angle. 

3 The local angle-of-attack in the r-s-t system 

a c is a measure of the force acting on the canard, 
expressed as an angle-of-attack. The inertial frame 
components of this vector are defined as axial, side, and 
normal canard angles-of-attack (viz., a cx , a cy , and a cz ). 

The “effective” canard angle-of-attack vector, a ce , 
is obtained by averaging a c through one period of the 
pitching sequence. The inertial frame components of 
this quantity are intended to describe the effect of a 
canard pitching sequence on the roll-averaged forces. 
Figure 41 shows results for a case for a = 0° and com- 
mand level 100%. The solid line is the C z obtained from 
the V2 simulation results (scale shown on the left verti- 
cal axis). The dashed line is the canard normal angle- 
of-attack, a cz , computed from the pitching sequence 
(scale shown on the right vertical axis). The dotted line 
is the canard pitching sequence (scale shown on the 
right vertical axis). The scales were chosen to show the 
excellent correlation between the V2 simulation results 
and the geometric quantity a cz (canard normal angle- 
of-attack). The side-force shows similar behavior. This 
indicates that the canard angle-of-attack vector is a use- 
ful estimator of forces on the canards. 

The effective canard angles-of-attack (i.e., roll- 
averaged inertial frame components of a c ) for the 
three canard pitching sequences at a = 0° are listed in 
Table 4. These numbers are computed from the pitch- 
ing sequence alone, and can be used as a first guess for 
trends in roll-averaged forces. An increase in command 
level gives mainly an increase in the effective normal 
canard angle-of-attack (a ), but also gives small 
increases in a cey and u cex . Increasing the command 
level from 100% to 200% increases the effective normal 
canard angle-of-attack by 27%. The corresponding 
increase for the computed normal force for the V2 solu- 
tion is also 27%. The corresponding numbers for effec- 
tive side canard angle-of-attack and V2 solution side 
force show increases of 19% and 26%, respectively. 
These are also in reasonably good agreement consider- 
ing that the side forces are small. 
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Figure 41. Canard NORMAL FORCE coefficient and angle-of-attack, a = 0°, c = 100% 


Table 4. Effective Canard Angles-of-attack* 



Command level 
0% 100% 200% 

a cex (axial) 
a cey (side) 
a cez (normal) 

2.7° 2.9° 3.3° 

0° 1.2° 1.4° 

0° 6.3° 8.0° 


Effective angles for three pitching sequences at 
free-stream angle-of-attack, a = 0°. 


In other words, lacking simulation (or experimen- 
tal) data for 200% command level, good estimates of the 
same are easily obtainable based on linear extrapolation 
as a function of effective canard angles-of-attack. Lin- 
ear extrapolation based on command level would fail, 
indicating an increase of 100% for both side force and 
normal forces. This illustrates that the roll-averaged 
forces are non-linear functions of the command level, 
and that effective canard angles-of-attack are better 
suited to aid interpolation and extrapolation of results 
between command levels. 


6. CONCLUSIONS 

A generic missile has been defined and analyzed to 
demonstrate the advantages and limitations of an overset 


grid approach to the aerodynamic performance predic- 
tion of rolling airframe missile systems. A computa- 
tional investigation has been carried out to determine the 
resolution required for accurate prediction of surface 
loads and to determine the relative aerodynamic signifi- 
cance of viscous and missile roll-rate effects. The aero- 
dynamic performance of the missile as a function of 
angle-of-attack and canard pitching sequence has also 
been considered. 

A moderate viscous grid system of 8 million points 
is a good compromise between solution accuracy and 
case throughput. Instantaneous and roll-average forces 
and moments from the moderate viscous grid are com- 
parable to those obtained in the high resolution baseline 
solution. Grid convergence of all forces and moments 
(except drag and rolling moment) were obtained. 

Significant viscous effects are apparent for the roll- 
ing airframe missile system considered. An inviscid 
solution has no capacity to account for interactions 
between the canard vortices and the missile boundary 
layer. The consequence of this inadequacy is incorrect 
positioning and strength of the inboard canard vortices 
and incorrect prediction of the secondary effects of the 
canards on fuselage and fin surface pressures. The 
inviscid results match well with the viscous solutions 
when the canards are in transition between maximum 
positive and negative deflection angles. During this 
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interval, the canard vortices are of minimal strength and 
the inviscid flow approximation is good. 

Variation of roll-rate in the range of 0 to 10 Hz does 
not have a significant affect on roll-averaged forces and 
moments for the FM-3 missile. The effects on instanta- 
neous loads are also minimal. The effects manifest as 
small roll-angle leads or lags in the load history. 

A single parameter method of approximation for 
changes in missile loads due to changes in canard pitch 
angle command level is presented. The approximation 
method helps interpret the physical significance of sim- 
ulation results and provides a powerful tool for interpo- 
lation of results between known data points. As such, 
the method is a quick and inexpensive means of supply- 
ing missing data, or for database expansion. 

An overset structured grid domain decomposition 
method enables accurate and efficient simulation of roll- 
ing airframe missile configurations that involve relative 
motion between system components. Case throughput 
rates are sufficient to contemplate aerodynamic database 
population with Navier- Stokes solutions. It is signifi- 
cant to note that time-accurate simulation for one rolling 
airframe case with OVERFLOW-D is much less expen- 
sive than generating the same data via corresponding 
quasi-static solutions for needed roll-angles and canard 
pitch positions. Each quasi-static solution costs a sig- 
nificant fraction of a single time-accurate rolling air- 
frame solution. A combination of time-accurate 
simulations for a relatively small number of pitching 
sequences and the canard angle-of-attack based approxi- 
mation method represent a powerful way to quickly 
populate a large aerodynamic data base. 
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